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Results of high resolution x-ray diffraction experiments are presented for single crystals of the spin 
gap compound BaCuSi2 0g in the temperature range from 16 to 300 K. The data show clear evidence 
of a transition from the room temperature tetragonal phase into an incommensurately modulated 
orthorhombic structure below ~100 K. This lattice modulation is characterized by a resolution 
limited wave vector q/c=(0,~0.13,0) and its 2 nd and 3 rd harmonics. The phase transition is first 
order and exhibits considerable hysteresis. This observation implies that the spin Hamiltonian 
representing the system is more complex than originally thought. 



INTRODUCTION 

BaCuSi206 is a quasi-2D compound composed of 
copper silicate layers separated by Ba ions. Within 
the Cu2Si40i2 layers, Cu 2+ ions are arranged in well- 
separated, vertical dimers 0, 0- Consistent with the 
dimerized structure, the material is observed to have a 
singlet ground state in zero magnetic field, with a large 
gap to the lowest excited triplet states . Magnetic fields 
in excess of H& ~ 23 T close the spin gap, such that cool- 
ing in a large applied field results in a state characterized 
by long-range magnetic order Q. At T — 0, a quan- 
tum critical point (QCP) at H c \ separates the quantum 
paramagnetic regime from the ordered state. Recent ex- 
periments probing the critical exponent associated with 
the approach of the phase boundary towards the QCP 
indicate that it is possible to describe the phase tran- 
sition in terms of Bose-Einstein Condensation (BEC) of 
delocalized triplets down to the lowest measured temper- 
atures of 30 mK [^, Q . Unfortunately, the large spin gap 
of this material currently precludes direct measurement 
of the magnetic structure in the ordered state. For this 
reason it is particularly important to have a detailed un- 
derstanding of the low-temperature crystal structure in 
zero magnetic field, since this will have important conse- 
quences for the nature of the high-field magnetic phase 
transition. 

To date, the crystal structure of BaCuSi20g has only 
been determined at room temperature and above. Sparta 
and Roth Q have found evidence for a subtle structural 
phase transition at 610 K, from a high temperature (HT) 
phase with IA/mmm (no. 139) symmetry, to a room 
temperature (RT) phase with IAi/acd (no. 142) sym- 
metry. However, recent heat capacity and susceptibility 
measurements indicate the presence of an additional first- 
order phase transition at approximately 100 K 0, which 
is also likely structural in origin. In this paper we present 
results of high resolution x-ray diffraction experiments 
performed at the Advanced Photon Source (APS) on sin- 
gle crystals of BaCuSi20g from 16 to 300 K. We find clear 



evidence for a first-order structural phase transition at 
approximately 100 K to a low temperature (LT) crystal 
structure, characterized by an orthorhombic distortion 
of the RT tetragonal structure, and the appearance of an 
additional incommensurate lattice modulation. We dis- 
cuss the origin of this effect and the consequences for the 
high field ordered state. 

EXPERIMENTAL METHODS 

Single crystal samples of BaCuSi206 were grown via 
a slow-cooling flux technique, as described elsewhere 
The crystals were well-formed, and had a small mosaic 
spread of ~ 0.02°. High resolution x-ray diffraction ex- 
periments were performed on the X-ray Operations and 
Research IBM bending magnet beamline y| at the APS, 
Argonne National Laboratory. This study was done with 
a focused beam of 20.016 keV x-rays. A Si(lll) analyzer 
crystal was used in order to suppress background and im- 
prove resolution. The sample was cooled using a closed- 
cycle refrigerator, with a base temperature of ~16K. A 
Si-diode sensor was placed about 1 cm below the sample 
to measure its temperature. 

Susceptibility measurements were made using a com- 
mercial Quantum Design MPMS5 SQUID magnetome- 
ter, for a field of 5000 Oe aligned parallel and perpen- 
dicular to the crystalline c-axis. Data were taken for 
both increasing and decreasing temperatures. Two dif- 
ferent instruments were used, with slight differences in 
the cooling rate and degree of undercooling on approach- 
ing a temperature set point. 

RESULTS 

The temperature dependence of the susceptibility of 
BaCuSi206 follows a form typical for weakly coupled spin 
dimer systems, and can be fit by a dimer model with a 
spin gap of 4.45 meV, as has been previously reported 
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FIG. 1: Temperature dependence of the susceptibility of 
BaCuSi2C>6 between 80 and 110 K, for an applied field of 
5000 Oe oriented parallel and perpendicular to the crystalline 
c axis (circles and triangles respectively). Data are shown for 
warming (solid symbols) and cooling cycles (open symbols), 
indicated by arrows. Jumps in susceptibility associated with 
the first order phase transition are indicated by vertical lines. 
Differences in the cooling rate result in the difference in the 
lower temperature jump for the two data sets. 



by Sassago and co-workers |3j • Superimposed on this be- 
havior, a weak jump in the susceptibility can be observed 
in the temperature range between 85 and 102 K (Figure 
1). The temperature at which the jump occurs depends 
on the thermal history of the sample, indicative of a first 
order phase transition in this temperature range. 

At room temperature, well above the transition ob- 
served in the susceptibility, charge Bragg peaks were 
observed consistent with the tetragonal structure pre- 
viously reported by Sparta 0- However, at 16 K the 
even-order (h,0,Q) Bragg peaks were observed to have 
split along the longitudinal direction as shown in Fig. 
2. In addition, incommensurate (IC) satellite peaks ap- 
peared on either side of the split Bragg peaks. A care- 
ful set of measurements revealed that the splitting is 
consistent with transformation twins of an orthorhom- 
bic (or weakly monoclinic) lattice. These twins are ro- 
tated with respect to each other about the 2-fold c-axis 
by ~ 89° as indicated in the inset to Fig. 2. From 
the splitting we found the degree of orthorhombicity de- 
fined as A = tr 2 — rr-s ~ 0.2%, where a D and b are 

the orthorhombic lattice parameters (ao < bo; Fig. 3). 
The c-axis lattice parameter also exhibits a weak reduc- 
tion (~0.5%) below the transition temperature. The 
incommensurate peaks are characterized by a wavevec- 
tor q IC = (0,- 0.13,0) and its 2 nd and 3 rd harmonics. 
These IC peaks are resolution limited in all three direc- 
tions implying a fully 3-dimensionally ordered modulated 
structure for the LT phase. The intensity of the funda- 
mental satellite is some ~ 10 5 times weaker than that of 
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FIG. 2: (Color on line). Reciprocal lattice scans (in red) along 
the b* axis showing the orthorhombic splitting of (8,0,0) 
Bragg peak at 16 K. IC diffraction peaks and their higher har- 
monics are observed on either side of the Bragg point. Note 
that the corresponding peaks on two sides are equidistant 
from (0,8,0), not from (8,0,0)t of the twin partner. Schematic 
diagram in the lower right corner illustrates split peaks from 
the twins. Scans along the orthogonal direction also reveal 
similar satellites (in blue; data were displaced down for clar- 
ity). In this case, the satellite peaks are associated with the 
(0,8,0)t twin. Central Bragg peaks have been scaled down by 
a factor of ~10 5 to display on the same scale. 



the (0,8,0) Bragg peak. 

The orthorhombic splitting and the intensity of the IC 
(0,8.13,0) peak were measured as a function of increas- 
ing temperature (Fig. 4). Both the splitting and the 
intensity of IC peak remain nearly constant up to 103.8 
K and disappear above that temperature. Within the 
resolution of our measurements we did not observe any 
change of q/c up to 103.8 K. An extensive search at 
110 K confirmed the absence of both the orthorhombic 
splitting and the IC modulation. On cooling, both the 
splitting and the IC peaks reappear below ~ 91K indicat- 
ing hysteretic effects, similar to what is observed in the 
susceptibility (Fig. 1). The hysteresis of the transition 
is spectacularly demonstrated by the behavior of the IC 
peak with thermal cycling through the transition (inset, 
Fig. 4). On warming from 16 K to 103.1 K, just below 
the transition, the peak is clearly observed (red circles). 
By raising the sample T to 104.6 K, just above the transi- 
tion, and subsequently cooling down to 100 K we were not 
able to recover the IC peak at 102.8 K (black squares). 
However, cooling down to 80 K and subsequently heating 
was sufficient to observe the peak up to 103 K on warm- 
ing again (blue triangles). From such measurements we 
deduced that the transiton is of first order, with transi- 
tion temperatures (T c ) of 104.1 ± 0.2K and 91K±1K on 
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FIG. 3: Temperature dependence of the lattice parameters 
for a warming cycle. 



warming and cooling, respectively. The observed transi- 
tion width is remarkably sharp (4£ ~ 0.003), which is 
a measure of high quality and purity of the bulk single 
crystal. 



DISCUSSION 

Low-T incommensurately modulated structures are 
not uncommon in silicates. Such modulations are usually 
due to a displacive phase transition of the lattice involv- 
ing rotations of SiC>4 tetrahedra without substantial in- 
ternal distortions of the tetrahedra themselves pS UlclllH . 
BaCuSi2 0e can be thought of as a 'vierer' single ring cy- 
closilicate 12], in which isolated S14O12 rings (composed 
of 4 corner sharing SiC>4 tetrahedra) are linked by CUO4 
groups in the afr-plane, and separated by Ba cations be- 
tween successive layers. It is likely that SiC>4 tetrahedra 
in this compound rotate and twist as a unit with respect 
to each other causing a lattice modulation involving both 
Ba and/or Cu atoms. The harmonic content (Figure 2) 
implies a subtle squaring-up of the modulated structure, 
but a detailed description of the atomic motions awaits 
a complete determination of the LT structure, which is 
beyond the scope of this initial survey. 

The susceptibility data shown in Figure 1 demonstrate 
that the change in crystal structure between the RT and 
LT phases of BaCuSi2C>6 affects the magnetism of the 
system, albeit very weakly. These changes are presum- 
ably associated with subtle changes in the superexchange 
parameters coupling the spins in the lattice. Without a 
detailed structural model, it is difficult to predict how the 
intradimer (J) and interdimer (J') exchange constants are 
affected, but we can anticipate that both will be modu- 
lated to some degree. We can therefore be fairly sure 
that the full spin Hamiltonian describing the system is 
slightly more complex than initially thought. We should 
note, however, that this structural study cannot provide 
an energy scale for the variations in the exchange par- 
tameters beyond suggesting that they are rather small, 
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FIG. 4: (Color on line.) Temperature dependence of the or- 
thorhombic splitting of the lattice (blue squares, left axis) and 
the intensity of fundamental IC modulation peak observed at 
(0,8.13,0) (red closed circles, right axis) measured on warm- 
ing. The intensity of the IC satellite on cooling is shown as 
open red circles. Note the hysteresis. Inset displays scans 
through the satellite region for three different thermal cycling 
profiles through the transition, discussed in the main text. 
The data in the inset correspond to the final temperature of 
the given sequence of temperature set points. 



given that the intensity of the incommensurate satellite 
peaks is several orders of magnitude lower than the Bragg 
peaks associated with the average crystal structure. 

The change in symmetry at the structural phase tran- 
sition can also affect the magnetism via spin-orbit cou- 
pling. The HT (IA/mmm) structure has a center of in- 
version symmetry between the two Cu ions that comprise 
a dimer unit. However, the RT (IAi/acd) and (presum- 
ably) LT structures do not, so additional terms D-Si x S2 
(where Si and S% label spins in one dimer unit) due to 
antisymmetric Dzyaloshinski-Moriya (DM) exchange are 
not forbidden in the spin Hamiltonian for these phases. 
This is the lowest order effect by which spin-orbit cou- 
pling can introduce terms to the spin Hamiltonian that 
break axial U(l) symmetry (an essential prerequisite 
for BEC |l4|'). since anisotropy in the (/-tensor simply 
rescales the critical field H c i = A/gfXB- In the RT struc- 
ture, copper dimers reside on axes of improper rotation 
4 (an inversion tetrad), which within a single CU2S14O12 
sheet closely resembles an axis of 4-fold symmetry. Ac- 
cording to the standard symmetry rules for DM exchange 
[loj l , the corresponding D- vector for the RT phase would 
then point along the 4-fold axis. Hence, for the RT struc- 
ture it is anticipated that the U(l) symmetry of the spin 
Hamiltonian will be preserved for magnetic fields aligned 
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parallel to the crystalline c-axis. Without a complete 
structural model for the LT phase we have less insight to 
the symmetry of the dimers. However, we note that the 
structural distortion is slight, which implies that there 
will be minimal changes in the D vector. 

In principle, the orthorhombicity of the LT structure 
can introduce an additional anisotropy to the spin Hamil- 
tonian via a second order effect associated with spin-orbit 
coupling. For a tetragonal structure, the x and y compo- 
nents of the interdimer exchange coupling J' are perforce 
equal, but for an orthorhombic structure this does not 
have to be the case. The resulting symmetric anisotropy 
is quadratic in the spin-orbit interaction, (i.e. only ap- 
pears in second order corrections) and is therefore even 
smaller than any DM terms. 

While the symmetry of the crystal lattice implies that 
the spin-orbit effects discribed above are possible, the 
magnitude of those terms must be measured by separate 
techniques. Ongoing ESR experiments indicate that the 
energy scale associated with the DM terms in both the 
LT and RT phases is negligible ^jf. Furthermore, the 
observed critical scaling exponents associated with the 
phase boundary imply that any U(l) symmetry break- 
ing terms that do enter the Hamiltonian must be on a 
significantly lower energy scale than the temperatures 
over which the magnetic ordering transition has been ob- 
served. 

Finally, we note that the observation of an incommen- 
surate structural modulation helps to understand results 
of preliminary high-field NMR experiments on this ma- 
terial. In these experiments, Stern and coworkers have 
studied the NMR spectrum of 29 Si nuclei at 0.04 K for 
fields just below H cl (23.4 T), and just above H cl (24.4 
Rather surprisingly, they find that the sharp 
NMR line below H c i is substantially broadened for fields 
above H c \. Within a simple interpretation of the pro- 
posed triplet BEC, one would instead anticipate a split- 
ting of the NMR lines corresponding to a doubling of the 
unit cell associated with the expected (it / a, it la) wave 
vector of the staggered magnetization[(|0]- The incom- 
mensurate structural modulation provides means to un- 
derstand this observation without necessarily having to 
introduce a more complex magnetic structure. At a tem- 
perature of 0.04 K for fields below H c \ there are essen- 
tially no thermally populated triplets, the uniform mag- 
netization is zero, and the Si nuclei all see the same local 
magnetic field resulting in a relatively narrow NMR line. 
However, for fields above H c \, there is a finite magneti- 
zation which causes a substantial Knight shift. Due to 
the incommensurate structure, each Si nucleus would see 
a slightly different magnetic field, broadening the NMR 
lincshape. That is to say that the broadened NMR lines 
seen for fields above H c \ do not necessarily correspond 
to an incomensurate magnetic structure, but more likely 
reflect the underlying incommensurate crystal structure 
of the LT phase. Further NMR experiments are called 



for to experimentally distinguish these two effects. 



CONCLUSIONS 

In summary, our experiments have unambiguously 
shown the presence of a low-temperature structural 
phase transition in BaCuSi2 0g, corresponding to an or- 
thorhombic distortion of the RT IA\/acd structure and 
accompanied by an incommensurate modulation. The 
transition is first order, having a substantial hysteresis 
spanning approximately 80 to 105 K, depending on de- 
tails of the warming or cooling cycles. The lower sym- 
metry and incommensurate modulation of the LT phase 
imply that the spin Hamiltonian describing the magnetic 
properties of the system is more complex than previously 
thought. ESR and inelastic neutron scattering experi- 
ments are currently in progress to determine the energy 
scale for additional relevant terms. Finally, the observa- 
tion of an incommensurate structural modulation at low 
temperatures provides important insight to recent high 
field NMR measurements for this compound. 
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